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The CALICE collaboration investigates different technology options for highly granular 
calorimeters for detectors at a future electron-positron collider. One of the devices 
constructed and tested by this collaboration is a lm3 prototype of a scintillator-steel 
sampling calorimeter for hadrons with analogue readout (AHCAL). The light from 7608 
small scintillator cells is detected with silicon photomultipliers. The AHCAL has been 
successfully operated during electron and hadron test-beam measurements at DESY, 
CERN, and Fermilab since 2005. One of main tasks for the successful data taking is 
to establish procedures for the equalization of cell responses and calibration of such a 
large number of channels. 

1 Introduction 
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The CALICE collaboration [T] is developing a hadronic calorimeter (HCAL) with very high 
granularity for future linear colliders (ILC, CLIC). The collaboration built a 1 m 3 physics 
prototype in 2005-6 [5] and currently CALICE is building an engineering prototype (3J- 

The AHCAL (Analogue Hadron CALorime- 
ter) 1 m3 physics prototype [2] with 7608 ac- 
tive readout channels has been in test beams 
at the CERN SPS in 2006 and 2007, fol- 
lowed in years 2008 and 2009 at Fermilab. 
The physics prototype is made of 38 lay- 
ers with scintillator tiles, interlaced with 16 
mm Fe absorber plates The dimension of 
one layer is 90 cm x 90 cm and it contains 
scintillator tiles with different granularity 
20 pes of 12x12 cm 2 tiles, 96 pes of 6x6 cm 2 
tiles and the first 30 layers have 100 pes of 
3x3 cm 2 tiles in the middle. The last 8 lay- 
ers have 6x6 cm 2 tiles in the middle instead. 
The active layers are referred to as modules, 
and the sum of active and passive material 
adds up to a total depth of 5.3 nuclear in- 
teraction lengths (A,). 

The active element of the readout is a 
compact photo-detector — a silicon pho- 
tomultiplier SiPM with 1156 pixels, each 

32/imx32/im in size. The AHCAL analogue board carries 1 ASIC for 18 SiPM channels. 
A second board contains control and configuration electronics and provides correct voltage 
for each SiPM. For the calibration a board called CMB is used, see Fig[T] 
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Figure 1: Sketch of one AHCAL module con- 
taining the scintillator tiles with SiPM read- 
out routed to the VFE electronics (on the right 
side. The CMB (on the left) provides UV LED 
light for calibration. 
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The AHCAL readout chain of this prototype contains scintillator tiles with embedded 
wavelength-shifting fibers and small SiPM photo-detectors. The electrical signal is adjusted 
by a preamplifier and a shaper and digitized by a 12-bit ADC. 

The variation of characteristics of the complete chain (gain, saturation) depends mainly 
on changes of temperature and operation voltage of the SiPM. For the correct offline re- 
construction of the energy deposition, the calibration runs have to be included into the 
data-taking process, since the condition inside the detector can change. For purposes of 
monitoring the long-term stability and performance of the photodctectors the calibration 
system with optical fibers for a UV LED light distribution was built [4]. A calibration and 
monitoring board (CMB) connected to each module distributes UV light from an LED to 
each tile via clear fibers. The LEDs are pulsed with 10 ns wide signals steerable in ampli- 
tude. By varying the voltage, the LED intensity covers the full dynamic range from zero 
to saturation (about 70 times the signal of a minimum- ionizing particles — MIP). Further- 
more, the LED system monitors variations of SiPM gain and signal response, both sensitive 
to temperature and voltage fluctuations. The LED light itself is monitored with a PIN 
photo-diode to correct for fluctuations in the LED light intensity. 



2 Calibration procedure 

For trustworthy data measurements with the AHCAL prototype we have to establish a 
reliable and robust calibration chain. This requires measurements with beam particles and 
with light from the LED monitoring system. The calibration chain is summarized in the 
following steps: 

• calibration of the cell response and cell-to-cell equalization; 

• monitoring of the SiPM gain and corrections for the non-linear response; 

• calibration to an energy scale (in GeV) with electromagnetic showers. 

For an energy deposit Ei in units of MIP of one single cell i (with registered signal Ai in 
ADC counts) we can write 

E * t MIP ] = A %n^ ■ /«t(^ tP«l), (1) 

where the first term is calibrated signal response in a given cell in case of a linear device. Due 
to the non-linearity in the SiPM response (finite number of SiPM pixels) , we need to add a 
correction function, which is described by the second term in the Eq. [T] This correction is 
a function of a number of SiPM pixels, A^[pix], firing for a single cell i. It is related to the 
ADC value for the cell, A,[ADC], and the corresponding SiPM gain, C,f ix [ADC]: 

Ai [pix] = Ai [ADC] /Cf x [ ADC]. (2) 

The first term of Eq. Q] expresses the equalization of all cell responses done by the cell- 
by-cell calibration, see the section |2~T1 The section [3] describes how to determine the gains 
for single cells with respect of temperature dependence. Methods how to correct for the 
non- linear SiPM response are introduced in the section 0] 
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2.1 Cell response equalization with MIPs 

The cell-by-cell calibration (equalization of all cell responses) is performed using minimum- 
ionizing particles (MIPs) provided by a broad muon beam, illuminating all cells in the 
detector. For each cell, a calibration factor, Cj MIP , is determined from the most probable 
value of the measured energy spectrum for muons in ADC units, which is extracted with a 
fit using a Landau function convoluted with a Gaussian, as can be seen for one cell in FigO 
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This fit accounts for the distribution of 
energy loss of muons in the scintillator tiles 
as well as for contributions from photon 
counting statistics and electronic noise. The 
combined systematic and statistical uncer- 
tainty for these fits was typically on the or- 
der of 2%. The muons are generally par- 
allel to the beam line and perpendicular to 
the detector front face. In this way all cells 
can be calibrated at the same time, mini- 
mizing the impact of temperature induced 
variations. 

The MIP calibration efficiency, i. e. 
the amount of successful MIP calibrations 

above the 0.5 x MIP threshold, is about 93% for the 2007 data set (CERN tests). The cali- 
bration fails for cells with broken electrical connection or very noisy SiPMs. The signal over 
noise ratio (N/S), defined as MIP amplitude over pedestal width (see Fig. [5]), is a measure 
for the separation of a MIP signal and noise. It is found to be ~10 in 2007. 



Figure 2: Single calorimeter cell response to 
muons with the corresponding fit (green line) 
and noise spectrum of the same cell. 



3 UV LED calibration and monitoring system 




Since the SiPMs are sensitive to changes 
in temperature and operation voltage, the 
LED system has been installed to moni- 
tor the stability of the readout chain in 
time. The calibration system needs suf- 
ficient flexibility to perform several differ- 
ent tasks. Gain calibration: we utilize the 
self-calibration properties of the SiPMs to 
achieve the calibration of an ADC in terms 
of pixels that is needed for non-linearity cor- 
rections and for a direct monitoring of the 
SiPM gain. Further we monitor all SiPMs 
during test beam operations with a fixed in- 
tensity light pulse. Saturation: we cross 
check the full SiPM response function by 
varying the light intensity from zero to the 

saturation level. Our calibration boards use an UV-LED as a light source for calibration. 
The UV-LEDs require a special driver [5] in order to make them shine fast (<10 ns) with 



Figure 3: Single photo-electron peak spec- 
trum taken with a SiPM, fitted with a multi- 
gaussian formula (red curve). 



3 



LCWS11 



■100, 



90 



S 80 



70 



CERN 07 
FNAL 08/09 



20 40 60 80 100 
measurement 



■100 



90 



'cS 80 



70 



CERN 07 
FNAL 08/09 



20 40 60 80 
measurement 



Figure 4: Gain calibration efficiency (left) and electronics inter-calibration efficiency (right), 
data taken at CERN in 2007 and at FNAL in 2008. More than 85.0 % of the channels could 
be monitored for gain and inter-calibration variation during these periods. 



an amplitude covering several orders of magnitude in light intensity. 
3.1 SiPM gain calibration 

For the extraction of the saturation curve / sa t (Eq. [1]) describing the SiPM non-linearity we 
need to know the SiPM gain for each cell (Eq. [5J . The gain is extracted from single photo- 
electron (p. e.) spectra taken in dedicated runs with low LED light intensity. LED light is 
necessary as the best determination of the gain requires a single photo-electron spectrum 
with a Poisson mean of about I .5 p.e. and the mean obtained from dark noise events is 
below 0.5 p.e. 

The SiPM gain, Gf, is the distance between two consecutive peaks in the single photo- 
electron spectrum. A typical gain spectrum is shown in Fig. |3l A multi-Gaussian fit is 
performed to the single photo-electron peaks to determine their average relative distance. 
The mean of each Gaussian function in the multi-Gaussian sum is left as a free parameter and 
the width is dominated by electronic noise, but for large number of pixels fired the statistical 
contribution becomes visible, which lead to an increase of the peak width. Accordingly, the 
width of each peak is left as a free parameter. The uncertainty on the gain determination 
is mainly due to the fit and is about 2 % for fits which pass quality criteria. 

SiPM gain measurements were repeated approximately every eight hours during test 
beam operation. The SiPM gain varies with temperature and has to be corrected for it. 
The efficiency of the gain extraction is defined as the number of successful fits in one gain 
run divided by the number of channels which can be calibrated. Figure [4] (left) shows the 
efficiency of the gain extraction for a series of runs taken at CERN and FNAL, respectively. 
Initial problems during the system commissioning phase led to low efficiency, but after 
commissioning a gain extraction efficiency of about 95 % per run has been achieved. The gain 
efficiency was also stable after transportation and throughout the FNAL runs. Combination 
of several gain runs yields calibration of more than 99 % of all cells. The remaining 1 % of 
cells are calibrated with the average of the module to which they belong. 
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3.2 ASIC mode inter-calibration 

The measurement of SiPM gain is performed with a special mode of the readout chip (high 
pre-amplification gain, CM mode). In contrast, physics data taking are performed with 
approximately ten times smaller electronic amplification (PM mode). The inter-calibration 
factor, li, of the chip gain between both modes is used to determine the overall SiPM 
calibration factor (used in Eq.EJ, C*f ix = Gfjh. 

The extraction of the inter-calibration coefficients depends on the linear response of the 
chip in both modes for an overlapping range of input signals. The input signal is provided 
by the LED system. The amplitude of the signal is varied within the linear range by varying 
the LED light intensity. The response in each readout mode is fit with a line, and the ratio 
between the two slopes is the inter-calibration coefficient for one given readout channel. 
Ideally, this factor should be a simple constant between the two chip readout modes, but 
it turns out to depend on the SiPM signal form due to the different shaping times in the 
two modes. For longer SiPM signals the inter-calibration is bigger than for shorter SiPM 
signals. The inter-calibration factors between the chip readout modes range between 4 and 
13. 

As with the gain, the inter-calibration extraction efficiency is influenced by the quality 
of the LED light distribution system. The inter-calibration coefficient extraction efficiency 
during the 2007 and the 2008 data taking periods are plotted in Figure |4] (right). After 
commissioning was completed, all channels with the exception of the 2 % inactive channels 
and the channels connected to a broken LED, could be inter-calibrated. For the missing 
inter-calibration values the average of the module to which a SiPM belongs is used instead. 

3.3 Temperature variation 

The SiPM gain and photo-detection efficiency arc temperature dependent. The product of 
the two determines the SiPM response, which typically decreases by 3.7%/K. A procedure 
has been developed to correct temperature-induced variations in the calorimeter response 
using temperature measurements in each module. This procedure and its stability will be 
described in more detail in [BJ. To account for the included temperature variations, the 
visible energy of each data set is scaled by -3.7%/K to the average temperature of the muon 
data used for calibration. 

4 SiPM non-linearity 

Due to the limited number of pixels and the finite pixel recovery time, the SiPM is an 
intrinsically non-linear device. The response function of a SiPM correlates the observed 
number of pixels fired, iV p i x , to the effective number of photo-electrons generated, N pc , 
including cross-talk and after-pulses. The response of a SiPM can be approximated by the 
function 

N pix = N tot -(l-e- N **/ N ™), (3) 

with iVtot the maximum number of fired pixels. This formula is a useful approximation 
for the case of uniform light distribution over the pixels and short light pulses. In the 
above approximation, one can extract a correction function / sa t(A p i x ) for the SiPM non- 
linear response as the residual to linearity of the inverted SiPM response function. For the 
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Figure 5: (Left) the SiPM non-linearity correction function, / sa t. The points are the tab- 
ulated data of iVpi x versus N pe for one SiPM in the AHCAL. (Right) Ratio of maximum 
number of fired pixels, N to t (mounted), measured with SiPM mounted on a tile to ./Vtot(barc) 
measured directly with bare SiPMs. 



analysis of measured data we use the tabulated values of a fit of the double exponential 
function, see Fig. [5] (left). 

For studies saturation behaviour of the 
individual SiPMs, they were not mounted 
on tiles, but were bare SiPMs. Therefore, 
all the pixels have been illuminated with 
light in a homogeneous way. The mea- 
surement results for all SiPMs installed in 
the AHCAL are given in [2]. The maxi- 
mum number of fired pixels (A to t(bare)) for 
each SiPM is extracted with a fit to the 
measured points using Eq. [3] The spread 
(RMS) in the values of N to t{bare) between 
all the curves is about 20%. SiPMs with 
N tot (bare) > 900 have been pre-selected. 
This ensures not too large variations in the 
non-linear response function of each device. 

Alternatively, A to t has also been ex- 
tracted using the AHCAL LED monitoring 
system from measurements with the SiPM 

mounted on a tile (N to t (mounted)). The saturation point in number of pixels is independent 
on the linearity of the light. Fig. [5] (right) shows the ratio of A^ot (mounted) to A to t(bare). 
The plot shows that the maximum number of pixels in the in-situ setup is on average 80.5 % 
of the value determined in the laboratory setup [7] with a wide distribution (RMS=9%). 
This factor is interpreted as geometric mismatch between the WLS fiber and the photode- 
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Figure 6: Relative uncertainty of the deter- 
mination of the saturation point for a single 
channel (RMS/N to t)- Temperature corrected 
data covers one full beam test period. 
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tector. The fiber has a 1mm diameter while the SiPM active surface area is lxl mm 2 ; the 
geometric ratio between areas is 79 %, in agreement with the measured value. Therefore, 
only a fraction of the SiPM surface is illuminated and the laboratory curves are re-scaled 
by the measured value of 80.5 % to correct for this effect before they are used to correct for 
the SiPM saturation. 

The uncertainty of the determination of 
the saturation point in-situ (N tot ) for a sin- 



gle channel is lower than 3%, if the LED 
light range properly covers the SiPM satu- 
ration region, and if this region is measured 
well below the ADC saturation. Unfortu- 
nately, these conditions are true only for a 
sub-sample of about 73% channels as can 
be seen in Fig. |6j 

The impact of the saturation correction 
is seen in Fig. [7] where the energy per hit is 
shown with and without the correction fac- 
tor f sat applied, for 30 GeV electromagnetic 




showers. Whereas the correction is negli- 
gible for low signal amplitudes, it becomes 

significant at larger amplitudes, resulting in Figure 7: Hit energy spectrum for 30 GeV 
a strong correction for the tail of the distri- positron showers in the AHCAL. Open circles 
bution. For the maximum energy deposited (dots) show the data before (after) correction 
in one cell for a 30 GeV the correction factor for the non-linear response of the SiPM. 
is / sa t(A 4 ) ~ 3. 



5 Validation of the AHCAL calibration 



The linearity of the calorimeter response for 
a large range of incident particle energies 
is a key feature, which allows for an im- 
portant test of the calibration chain. Elec- 
tromagnetic showers offer the most rigorous 
test for non-linearity correction, since the 
energy deposited per single tile in an elec- 
tromagnetic shower is larger than that in a 
hadronic shower for the same particle en- 
ergy. Figure [8] shows the hit energy spec- 
trum of a 40 GeV positron shower compared 
to the spectra of 40 GeV and 80 GeV pion 
showers. The positron shower clearly has 
more hits with high energy deposition, even 
when the total particle energy is only half 
that of the pion. 

More results about the analysis of the 
electromagnetic showers measured by the AHCAL can be found in [8]. 
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Figure 8: Hit energy spectrum for 40 GeV 
positron showers compared to that of 40 GeV 
and 80 GeV pion showers. 
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